Quasiparticle Charge" in Superconductors: Effect of Mott Physics or Hidden Order 
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The renormalization factor, dubbed "quasiparticle charge" Ze, of the coupling between the su- 
percurrent and quasiparticles was analyzed in the context of high temperature cuprates in a recent 
paper by loffe and Millis.|l| They observed that Ze in cuprates deviates from the BCS value [Z^ — 1), 
which was interpreted as the proximity effect near a Mott insulator. Here we show that the de- 
viation from Ze = 1 can occur, in general, even in the absence of quasiparticle interactions, when 
the superconducting order coexists with another order parameter with the same internal symmetry. 
As an example, we compute the coefflcient of the linear temperature dependence of the superfluid 
density when the d-wave superconducting state coexists with the orbital antiferromagnetic state 
(d-density wave), and find that Ze varies from 1/^2 to 1. 

PACS numbers: 74.20.-z, 74.25.-q 



I. INTRODUCTION 

The behavior of charge transport in both the normal 
and superconducting phases of cuprates is one of the puz- 
zling issues in high temperature superconductors. In a 
recent paper, lofFe and Millis discussed the effect of prox- 
imity to a Mott insulator on charge transport in the su- 
perconducting state of cuprates. They pointed out 
that the renormalization factor, dubbed "quasiparticle 
charge" Z^,, of the the coupling between the supercur- 
rent and quasiparticles is renormalized by quasiparticle 
interactions and contains useful information about strong 
correlations near a Mott insulator. || 

Even though the real quasiparticle charge is not con- 
served in superconductors and is not the conven- 
tionally defined charge, we will still use "quasiparti- 
cle charge" to call the renormalization factor mentioned 
above for convenience. As pointed out in the past 
10, |2| ^, Ze appears in the linear temperature depen- 
dence of the superfluid density in d-wave superconduc- 
tors (dSC). The low temperature superfluid density of 
d-wave superconductors is given by 
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Here vp is the Fermi velocity and va the gap velocity 
defined as va = Aoa/-\/2, where Aq is the maximum 
value of superconducting order parameter. 

The behavior of the superfluid density in cuprates has 
been an important subject of debate. Lee and Wen ob- 
served that, in cuprates, a does not strongly depend on 
the doping concentration, x, while the zero temperature 
superfluid density is proportional to .t.|^, ^ This would 
be consistent with the phenomenology that the quasipar- 
ticles in the superconducting state of cuprates still cou- 
ples to external electromagnetic field in the usual BCS 



manner even in the presence of strong correlation near a 
Mott insulator. Q On the other hand, two different strong 
coupling approaches to the t-J model lead to different 
behaviors of the superfluid density. The U(l) mean field 
theory predicts a oc as well as Ps(0) oc x imply- 
ing that the "quasiparticle" charge goes to zero as the 
Mott insulator is approached while the SU(2) gauge the- 
ory lead to a ^ 0{1) with the same ps(0).^, ^ 

Later on, Millis et al pointed out that a is in 
general renormalized by quasiparticle interactions even 
in the conventional d-wave superconductors. Millis et al 
further argued that the deviation from the standard BCS 
value Zf, = 1 should be regarded as an evidence of the 
strong correlation present near a Mott insulator. Re- 
cently, loffe and Millis extracted the values vp and 
VA from the angle resolved photoemission spectroscopy 
(ARPES) and the specific heat data of cuprates, and ob- 
tained the value of Z^. for several samples of cuprates. 
They found that Z^ varies from 0.6 to 1 in cuprates, and 
argued that the deviation from the BCS value {Z^ = 1) 
is due to strong quasiparticle interactions near a Mott 
insulator. 

In this paper, we show that the deviation from the 
BCS value of Zg can also occur within the weak cou- 
pling approach without taking into account the quasipar- 
ticle interactions, when a superconducting state coexists 
with another ordered state with the same internal sym- 
metry. In this case, the value of Z^ contains informations 
about the additional order parameter. As an example, 
we consider a system where the d-wave superconducting 
order coexists with the commensurate orbital antiferro- 
magnetic order (d-density wave).[Q The quasiparticles in 
this state have the Dirac spectrum near the nodes as the 
case of dSC and the superfluid density has the linear 
temperature dependence. By computing the coefficient 
of the linear temperature dependence of the superfluid 
density, we find that Z^ in Eq]| varies from l/^/2 to 1. 
The deviation from the BCS value {Z^ = 1) in this case 
is not due to the quasiparticle interaction, but due to the 
existence of the additional order parameter. The super- 
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fluid density of the dSC coexisting with the dDW order 
was previously studied |l^ , but the question of the 
renormalization factor Ze was not addressed. Although 
we studied a specific system, our results can be applied 
to more general situation where a superconducting or- 
der coexists with another order parameter with the same 
internal symmetry. 

In the next section, we compute the current-current 
correlation function in the coexisting dSC and dDW 
states. Those who are not interested in the details can 
skip the next section and jump to section III. In the third 
section, we provide the analytical expression of a and ex- 
tract the value of Zg. We compare our results with the 
consequences of other existing theories and discuss fur- 
ther implications in the final section. 



II. SUPERFLUID DENSITY IN THE 
COEXISTING DDW AND DSC STATES 

In the phase where the dSC and dDW coexist, the 
quasiparticles keep the Dirac spectrum as the case of 
dSC. In particular, when the chemical potential /i is zero 
or /Lt < r <C min(Ao,Wo), the low energy nodal spec- 
trum is well described by changing the gap velocity va 
to WA = i/i^A + where vw = Woa/^/2 and Wq is the 
maximum value of the dDW order parameter. In this 
case, the gap velocity measured in ARPES will corre- 
spond to Va- The density of states will be also modified, 
but in a trivial fashion; va is replaced by va- Thus one 
may expect that the coefficient of the linear temperature 
dependence of the superfluid density is simply given by 
a = (ln2)vp/ {2ttva)- If this is the case, we will not 
be able to extract an independent information about vw 
and practically = 1; experiments will not be able to 
tell the difference between the simple dSC state and the 
coexisting dSC/dDW states. We show below that this is 
not the case. Z^ has an independent information about 
Vw and in general it is not unity. The deviation from 
Ze = 1 is in fact due to the existence of the dDW order 
parameter. 

We derive the result described above by directly 
computing the current-current correlation function in a 
mean-field theory. We start from the mean-field Hamil- 
tonian given by 

H ^(ek-/i)4,^Ck,^-f ^iWk4_^Ck+Q,^ 
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where 
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Here Ak and iWk are the order parameters of the dSC 
and dDW respectively. Using the Nambu basis defined 
as follows, 
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the current operator can be written as 

ja(q) = i^sin(fc„-hg„/2)*J^(p3/)^'k+q 



iWo sin {ka + ga/2)1'J,(ip2cr3)«'k+q, (8) 



where a — x,y. Here p's and cr's are the Pauli matrices 
and piaj = pi® aj. 

The superfiuid density can be computed from the 
current-current correlation function. llql 



p,iT) = i-K,) - lim lim A,,(q, i^), (9) 

where {—Kx) is the diamagnetic contribution and 
AxK (q, v) is the paramagnetic current-current correlation 
function. The diamagnetic part can be computed as 

{-Kx) = 2_, — Q smfc^ 

k 

X ( tanh ^ + tanh :^\lO) 

where 



Eik,2k = -^(l^fcTAi)^ + A2, 
^k = Jel + Wl 



(11) 



The paramagnetic contribution contains the information about and can be obtained from 

Axxiq., ii^ni) = t'^T^^sinkx sin (kx + qx)Tr[G{k,iujn){p3l)G{'k + q,iu!„ -f ij/™)(p3/)] 

71 k 

-fW^o^T^ ^ sinfca; sin (fca; + qx)Tr[G{'k,iujn){p2<J3)G{'k + q,iu!n + mn){p2<J3)], (12) 
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Here Vm = 2rmTT and ujn = (2n + 1)ttT are the bosonic and fermionic Matsubara frequencies. Notice that the 
cross-terms proportional to tWo cancel out each other. 

After the frequency sum, the correlation function has the following form. The first term, A^^, of EqjlJ which comes 
from the conventional current proportional to t is obtained as follows by taking v first and q ^ later. 
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where f{x) = l/(e="/^ + 1) is the For mi function and 
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On the other hand, the second term, Aj^^, of Eq.[l^ 
comes from the additional contribution to the current 
due to the existence of the dDW order parameter. After 
taking — > first and q ^ later, we get 
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where 
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Therefore, the supcrfluid density can be obtained from 



Eq.|, where A^^, = A^^ + A^ 
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III. VALUE OF THE "QUASIPARTICLE 
CHARGE" FROM SUPERFLUID DENSITY 

At low temperatures, the leading paramagnetic contri- 
bution comes from the nodal quasiparticles and we can 
obtain the low energy quasiparticle dispersion by expand- 
ing tk near /i, and Afc, Wk near the node as follows. 



Wk = vwv-, 



(18) 



where p+ = (p2,-|-pj,)/\/2 andp_ = (pa;—pj,) /a/2 with the 
momentum, p, measured from k^?. Here, vp = 



VA — Aoa/v^, and vw = Woa/V2. When ^ = or 
fi < T, the quasiparticle dispersion can be well described 

by 



^lk,2k 



2 2 
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Notice that the corresponding density of states has the 
following form. 



N{E) = 



E 



2nvF\/vA + ^' 



(20) 



Therefore, the specific heat as well as the measured gap 
velocity in ARPES will be given by the combination of 
Vw and ua through y/v'^ + Vy^. Now it is clear that 
the leading linear temperature dependence of the super- 
fluid density comes from the terms with df{Eik)/d{Eik) 
factor. The coefficient of the linear temperature depen- 
dence, a, is obtained as 



In 2 
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(21) 
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Comparing Eq.|l| and Eq.p|, an d replacing va by the 

we find 



measured gap velocity, 



2{vl 



(22) 



where we neglect 0{W'^/t^) terms. From the above re- 
sult, one can easily see that 1/2 < < 1. If Wq ^ Aq, 
then 1/2, while in the opposite case Aq ^ Wq 

or Wq ^ 0, we get ^ 1. Therefore, varies from 
1/V2{= 0.71) to 1 which is curiously close to the values 
in cuprates obtained by loffe and Millis 

In the discussion above, we assumed /i < T. In the 
opposite limit, /i > T, the low energy quasiparticle dis- 
persion is better described by 



Efk = {vFP+f + [vAP-f 



(23) 



which does not have any information about Wk , the dDW 
order parameter. This is due to the fact that the quasi- 
particles recognize the presence of the dDW gap only 
when their energy scale becomes larger than the chemi- 
cal potential. Carrying out the explicit computation of 
the current-current correlation function, we verified that 
the coefficient of the linear temperature dependence of 
the superfluid density is given by Eq.| with = l.(l4| 
Therefore, there is a cross-over in the value of Z^, to unity 
when ^ > T while Z^ ^ I ior p < T. At the half-filling 
(/i = 0), Ze is in general not unity. At any rate, if is 
smaller than T in the experimentally relevant tempera- 
ture range, the value of Ze is given by Eq.|2^. 



IV. CONCLUSION 

In the strong coupling regime, the coupling between 
the supercurrent and the quasiparticles can be strongly 



renormalized by quasiparticle interactions and lead to the 
deviation of Z^ from the BCS value {Zf. = 1). Here we 
show that the deviation from Z^, — I can also occur in 
the weak coupling regime when the superconducting or- 
der coexists with an additional order parameter with the 
same internal symmetry. We have shown that in the case 
of the coexisting dSC and dDW states, the "quasiparti- 

cle charge" is given by Ze = f 1 - 2{v'i+v-i,) ) "^^^^ 
the chemical potential ^ is zero or fi < T, where vw and 
VA are the gap velocities determined by the maximum 
gaps of the dDW and dSC order parameters. Thus Z^ 
varies from l/-\/2 to 1. The deviation of Z^ from the 
unity is a good measure of the additional coexisting or- 
der parameter. Our result of Ze ^ 1 is also applicable 
to the systems where a superconducting order parameter 
coexists with another order parameter with the same in- 
ternal symmetry; e.g., p-wave superconductor coexisting 
with p-density wave, although the expression of Ze would 
be different for different systems. 



Our results suggest that the determination of the value 
of Ze alone will not sharply distinguish two possibilities; 
Mott physics and the existence of a hidden or an addi- 
tional order parameter. However, if there exists a su- 
perconductor far away from a Mott insulator while it 
exhibits Ze < 1, one would strongly suspect that there 
might exist a hidden order parameter. The direct rel- 
evance of our results for the coexisting dSC and dDW 
states to cuprates has to discussed in conjunction with 
other experimental data and it is beyond the scope of 
this paper. 
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